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Outgoing  longwave  radiation  (OLR)  simulated  by  the  GLAS 
general  circulation  model  is  validated  with  that  derived 
from  measurements  made  by  polar  orbiting  satellites.  The 
mean  observed  and  simulated  OLR  fields  are  compared  for  spa¬ 
tial  and  seasonal  variations  by  expanding  them  into  spheri¬ 
cal  surface  harmonics.  The  student’s  *t’  test  is  used  to 
determine  the  significance  of  differences  between  the  simu¬ 
lated  and  observed  fields.  .rr?.  '//rr 

The  global  mean  and  standard  deviation  of  the  simulated 
geographical  fields  are  respectively  about  20  Wm~*^smaller 
and  10  Wm~*  larger  than  the  corresponding  values  of  the  ob¬ 
served  field.  The  major  fraction  of  the  variability  in  the 
simulated  field  is  due  to  a  sharp  meridional  gradient  in  the 
field.  The  seasonal  variation  of  the  global  mean  and  stan¬ 
dard  deviation  of  the  simulated  fields  are  in  good  agreement 
with  observations.  Correlation  coefficients  between  the  ob¬ 
served  and  simulated  fields  as  a  function  of  spatial  scales 
show  that  the  phase  relationship  for  large  spatial  scales  is 


,very  good  -for  January  and  July  Months  but  -fair  for  April  and 


October.  In  the  tropics,  the  differences  between  sieulated 
and  observed  OUR  Means  are  not  significantly  different  f roe 
zero,  except  over  deep  convection  regions  (Asiatic  aonsoon, 
Anazon,  Central  Africa)  where  the  Model  convective  clouds  do 
not  interact  with  radiation.  The  differences  in  the  Middle 
and  high  latitudes  are  highly  significant.  More  so  in  the 
southern  heel sphere  than  in  the  northern  heei sphere. 
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Validation  of  any  general  circulation  modal  (GCM)  out¬ 


put  involve 


■form  of  comparison  of  simulated  data  ver¬ 


sus  observed  data.  In  some  sense  the  more  accurate  and 
dense  our  observations  are,  the  more  complete  our  validation 


can  be.  However ,  our  knowledge  o-f  such  basic  meteorological 
parameters  as  temperature,  pressure,  and  winds  is  limited 
due  to  sparse  data,  particularly  over  the  oceans.  Compound¬ 
ing  the  validation  problem  is  the  accuracy  of  these  data, 
which  may  contain  some  errors  as  large  as  lOX. 

Fortunately  the  measurement  of  at  least  one  atmospheric 
parameter  by  our  present  observation  network  is  both  uni¬ 
formly  dense  and  sufficiently  accurate  for  our  purposes. 
Twice  daily  measurements  of  outgoing  longwave  radiation 
(OLR)  for  the  past  10  years  (June  1974  to  May  1984,  except 
March— December  1978)  are  available  on  a  2.5  X  2.5  degree 
global  grid.  The  measurements  were  made  by  NOAA  satel¬ 
lites.  Specific  information  on  the  observed  OLR  data  is 
given  in  Section  2. 

Validation  of  any  model's  OLR  is  r<  ally  a  validation  of 
how  well  the  model  simulates  the  quantities  which  govern 
OLR.  These  include  the  temperature  of  the  earth-atmosphere 
system  (inherent  in  this  is  the  temperature  of  the  radiating 

surface,  be  it  land,  sea,  ice,  or  snow),  the  distribution  of 
cloudiness  and  the  distribution  of  the  major  OLR  absorbing 


and  emitting  gases  (water  vapor,  ozone,  and  carbon  diox¬ 
ide).  The  eost  important  quantities  -for  determining  OLR  are 
the  distribution  and  type  o-f  clouds.  As  a  simple  example, 
an  outbreak  of  cold  air  over  the  central  plains  of  the  U.S. 
would  cause  a  slight  decrease  in  the  OLR  from  that  region, 
while  a  cloud  layer  would  cause  a  large  reduction  in  OLR. 
To  date,  validation  of  the  OLR  -field  has  been  unsatisfactory 
-for  any  general  circulation  model  because  virtually  all 
6Cn*s  reproduce  the  critically  important  cloud  type,  amount, 
and  distribution  so  poorly.  However,  we  -feel  rapid  progress 
in  improvements  of  6Cn*s  in  this  area  justifies  OLR  valida¬ 
tion. 

It  seems  apparent,  then,  that  a  model  incorporating 
some  type  of  cloud  interaction  capability  will  yield  better 
results  when  using  OLR  as  a  parameter  for  validation,  than  a 
model  without  the  capability.  For  this  reason  we’ve  chosen 
the  GLAS  (Goddard  Laboratory  for  Atmospheric  Sciences)  gen¬ 
eral  circulation  model  for  our  study.  The  model’s  treatment 
of  cloud-radiation  interactions  is  discussed  in  detail  in 
Section  3. 

The  earliest  attempts  at  validation  of  the  output  of 
GCH’s,  based  on  zonally  averaged  data  fields  and  vertical 
cross-sections  of  the  atmosphere,  involved  rather  qualita¬ 
tive  and  subjective  comparisons.  While  these  comparisons 
are  a  useful  starting  point,  a  rigorous  scientific  valida¬ 
tion  requires  more  quantitative  and  objective  procedures. 

The  purpose  of  this  study  is  to  validate  with  satellite 


observations  th»  seasonal  variations  of  OLR  sieulated  by  the 
6LAS  model.  The  data  sets  are  composed  of  mean  annual  vari¬ 
ations  calculated  from  9  years  of  satellite  observed  OLR  and 
2  years  of  OLR  simulated  by  the  model.  We  represent  the  two 
data  sets  as  a  finite  sum  of  spherical  surface  harmonics, 
somewhat  along  the  lines  of  Christidis  and  Spar,  (1981). 
This  decomposition  allows  a  comparison  between  the  two 
fields  of  the  annual  variations  of  the  amplitudes  of  the 
various  components  and  their  correlations.  A  discussion  of 
our  application  of  spherical  harmonics  is  included  in  sec¬ 
tion  4.  The  6LAS  model’s  ability  to  simulate  OLR  is  evalu¬ 
ated  in  Section  5  with  a  variety  of  comparisons  of  the  sim¬ 
ulated  vs  observed  OLR.  These  include  a  season-by-season 
subjective  comparison  of  the  two  OLR  fields,  as  well  as  a 
brief  analysis  of  fields  representing  differences  between 
the  two.  Also  included  are  a  comparison  of  the  average 
monthly  means  and  a  plot  of  values  calculated  using  a  stu¬ 
dent’s  t-statistic. 


2. 1  Background 

NOAA  polar  orbiting  satellites  have  provided  measure- 
eents  of  outgoing  longwave  radiation  (OLR)  for  approximately 
9  years.  We  now  examine  the  consistency  and  accuracy  of 
these  data  to  establish  their  usefulness  in  GCM  validation 
studies. 

The  period  of  observations  spans  9  1/2  years,  fro*  June 
1974  through  November  1983,  with  a  10  month  gap  during  March 
-  December  1978.  The  remaining  nearly  9  years  of  data  have 
been  gathered  by  seven  different  satellites.  The  first  four 
satellites,  NOAA  2,3,4,  and  S  are  members  of  the  NOAA  scan¬ 
ning  radiometer  (NOAA  SR)  series,  and  they  provided  the 
largest  continuous  portion  of  the  record  (Gruber  and  Win¬ 
ston,  1978).  The  10  month  gap  began  with  the  failure  of 
NOAA  5  and  ended  with  the  launch  in  January  1979  of  the  Ti¬ 
ros  N  satellite.  The  Tiros  N  was  replaced  in  February  1980 
by  NOAA  6  which  was  succeeded  by  NOAA  7  in  September  1981. 
These  satellites  incorporated  two  important  types  of  changes 
that  affect  the  homogeneity  of  our  datai  1)  different  equa¬ 
tor  crossing  times,  and  2)  different  water  vapor  window 
channels  used  for  measurement  of  OLR. 

2.2  Equator  Crossing  Times 

Table  1,  extracted  from  a  table  by  Gruber  and  Krueger 


(1984),  shows  ths  various  satsllitss  and  thsir  associated 


equator  crossing  times  givsn  in  local  sun  time. 


NOAA  SR  Sari  as  9*00  am  and  pm 
Tiros  N  3*30  am  and  pm 
NOAA  6  7*30  am  and  pm 
NOAA  7  2*30  am  and  pm 


Sampling  at  different  timms  of  thm  day  introduces  a  de- 
gram  of  inhomogmnsi ty  into  thm  data  record.  Gruber  and 
Krueger,  (1984)  showed  that  the  12-hour  global  average  OLR 
difference  varies  from  about  1-2  Wm~a  for  sunrise/sunset 
measurements  to  about  5-6  Wm~a  for  mid-afternoon/night  meas¬ 
urements.  From  this  they  infer  a  diurnal  global  average  OLR 
cycle  having  an  afternoon  maximum  and  an  early  morning  mini¬ 
mum  with  an  amplitude  of  about  6-7  Wm~a. 

To  minimize  the  effect  of  the  diurnal  variation  on  the 
record  the  average  of  the  two  daily  measurements  is  used.  A 
comparison  of  the  zonal  profiles  of  averaged  annual  OLR 
measured  by  each  satellite  shows  only  small  differences 
(Gruber  and  Krueger,  1984).  Gruber  and  Krueger  point  out 

that  these  differences,  on  the  order  of  a  few  Hm~a,  may  be 
attributable  in  large  part  to  inter-annual  variability  of 


th#  QLR,  and  do  not  detract  significantly  fro*  the  assump¬ 
tion  of  homogeneity  in  the  data. 

2.3  Procaduras 

Tha  original  NOAA  8A  sari as  window  channel  spanned  the 
region  froa  10.3  to  12.3  Mb.  On  Tiros  N  and  NOAA  &  the  win¬ 
dow  channel  was  narrowed  to  10.3  to  11.3  mbv  while  on  NOAA  7 
the  infrared  channel  used  for  OLA  measurements  covered  11.3 
to  12.3  Me  (Janowiak  at  al..  1903). 

In  addition  to  variations  In  the  window  channel,  anoth¬ 
er  factor  that  went  through  several  modifications  was  the 
algorithm  used  to  derive  total  OLA  from  observations  in  a 
narrow  band.  Hark  at  al . .  (1962)  found  the  relationship  be¬ 
tween  the  narrow  band  measurement  and  total  outgoing  long¬ 
wave  emittance  to  be  linear  in  a  sample  of  99  different  at¬ 
mospheres.  This  linear  relationship  was  used  on  data  gath¬ 
ered  by  the  NOAA  8R  series  of  satellites.  Then,  Abel  and 
Gruber  (1979)  showed  that  a  non-linear  algorithm  provided  a 
better  fit.  Their  relation  was  used  to  convert  the  entire 
record  then  existent.  Versions  of  this  new  non-linear  rela¬ 
tion  were  used  for  the  data  collected  by  the  different  win¬ 
dow  channels  of  the  Tiros  N  and  NOAA  6  satellites. 

However,  Ellingson  and  Ferraro  (1983)  in  comparing 
fluxes  from  the  NOAA  satellites  with  observations  from  the 
NXHBUB  6  Earth  Radiation  Budget  experiment  found  through 
their  radiance  calculations  that,  except  for  clouds  below 
730  mb,  the  relationship  derived  by  Abel  and  Gruber  over- 


NtiMtfd  the  OLR  -from  partial  cloud  covtr.  Th*  Magnitude 
of  the  errors  varied  Mith  cloud  type  and  aeount  but  in  com 
cases  Mas  as  large  as  2%  of  the  planetary  outgoing  flux. 
The  subsequent  NOAA  7  data  used  an  updated  algorithm  based 
on  the  work  of  Ellingson  and  Ferraro. 

Finally,  Ohring  et  al .  (1984)  reviewed  the  procedure 
for  converting  window  measurements  to  total  OLR  and  found 
that  the  record  based  on  the  previous  methods  contained  a 
positive  bias  of  about  13  Wm~a.  Their  work . involved  an 
analysis  of  simultaneous  measurements  from  the  NIMBUS  7  sat¬ 
ellite  of  the  total  longwave  emittance  and  the  radiance  in 
the  10—12  Mn  window  measured  by  the  Thermal  Infrared  Radio¬ 
meter  (THIR) .  The  new  algorithm  developed  by  Gruber  and 
Krueger  reduced  the  global  average  OLR  by  approximately  11 
Nm~3  (Gruber  and  Krueger,  1984).  By  relating  the  THIR  win¬ 
dow  radiance  to  the  window  radiances  of  each  of  the  satel¬ 
lites  participating  in  this  data  collection,  Gruber  and 
Krueger,  (1984)  also  derived  appropriate  algorithms  to  en¬ 
able  a  final  update  of  the  entire  record. 

Because  of  the  extensive  work  and  multiple  refinements 
in  the  algorithms  that  calculate  total  OLR  from  window  meas¬ 
urements  described  above,  we  believe  the  data  set  to  be  a 
good  representation  of  the  actual  total  OLR  field. 


evolved  from  the  3-livtl  GCM  developed  at  UCLA  by  A.  Arakawa 
and  Y.  Mintz  (Arakawa,  1972). 

The  GLAS  GCM  is  a  9-level  primitive  equation  model  in 
sigma  coordinates.  Each  of  the  nine  vertical  levels  has  the 
same  sigma  thickness.  The  total  mass  of  the  model  atmos¬ 
phere,  which  in  pressure  units  is  surface  pressure  minus 
pressure  at  the  top,  is  divided  into  nine  equal  layers.  The 
upper  boundary  of  layer  1  is  at  10  mb.  A  vertical  differ¬ 
encing  scheme  is  employed  based  on  the  work  of  Arakawa  and 
Lamb  (1977)  and  Arakawa  and  Suarez  (1983).  Horizontal  grid 
spacing  is  four  degrees  of  latitude  by  five  degrees  of  lon¬ 
gitude.  Variables  are  horizontally  staggered  according  to 
the  B-grid  scheme  of  Arakawa  and  Lamb  (1977).  To  maintain 
computational  stability,  selected  fields  are  filtered  in  the 
longitudinal  direction  near  the  poles  with  a  Fourier  filter 
based  also  on  the  design  of  Arakawa  and  Lamb  (1977).  The 
Fourier  filter  is  applied  to  the  zonal  wind,  which  is  used 
to  compute  the  zonal  mass  flux,  and  to  the  zonal  component 


), 


‘•W 


o-f  the  pressure  gradient  -force 


A  recently  modified  planetary  boundary  layer  (P6L)  pa¬ 
rameterization  (Randall,  1982a)  has  been  included  which  in¬ 
corporates  an  improved  surface  -flux  parameterization.  How¬ 
ever,  even  with  these  improvements  the  model  simulates  plan¬ 
etary  albedo  as  0.40  compared  with  an  observed  planetary  al¬ 
bedo  o-f  approximately  0.30  (Shukla  et  al . .  1981). 

The  GLAS  model  contains  sophisticated  routines  -for 
dealing  with  radiation  in  both  clear  and  cloudy  skies  over 
the  entire  spectral  range.  However,  since  this  study  -focus¬ 
es  on  the  outgoing  longwave  radiation,  only  the  model's 
treatment  of  longwave  radiation  will  be  addressed. 

3.2  Longwave  Radiation 

The  most  important  absorbers  o-f  longwave  radiation  in 
the  atmosphere  are  water  vapor,  carbon  dioxide,  and  ozone. 
The  GLAS  model  uses  the  Wu-Kapl an  radiation  routine  (Krish— 
namurthy,  1982),  which  is  summarized  here,  in  its  treatment 
o-f  these  absorbers. 

3. 2. 1  Absorbers 

The  absorption  o-f  radiation  by  atmospheric  gases  is 
spectrally  dependent.  The  spectral  region  where  the  absorp¬ 
tion  occurs  is  determined  by  the  molecular  structure  o-f  the 
gases.  Radiation  is  absorbed  in  spectral  bands  which  con¬ 
sist  o-f  many  individual  narrow  lines.  Water  vapor  has  a 
strong  vibration-rotation  band  at  6.3  Mm  and  pure  rotation 


.Sm *.  J 
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bands  at  wavelengths  greater  than  20  Mn.  In  the  Wu-Kaplan 
radiation  routine,  the  temperature  dependence  of  absorption 
coefficients  iti  the  6.3  Mn  water  vapor  vibration— rotation 
band  is  not  taken  into  account,  as  the  radiative  transfer  in 
that  band  is  relatively  small.  The  most  important  carbon 
dioxide  band  is  a  vibration-rotation  band  at  14.7  Mr,  while 
ozone  absorbs  intensely  near  9.6  Mn  with  other  bands  at  4.7 
Mn  and  14.1  Mn.  In  the  regions  between  the  molecular  ab¬ 
sorption  bands,  atmospheric  gases  are  relatively  transpar — 
ent.  However,  water  vapor  has  an  absorption  continuum  in 
the  atmospheric  window  region  between  approximately  8  Mn  and 
13  Mm. 


3.2.2  Transmission  Functions 

Of  the  three  major  absorbers,  detailed  calculations  are 
carried  out  for  water  vapor  only,  while  pre-calculated  quan¬ 
tities  such  as  transmission  functions  are  used  for  carbon 
dioxide  and  ozone.  Because  of  the  coniplicated  structure  of 
the  spectral  absorbing  bands,  it  is  a  formidable  task  to 
carry  out  a  line-by-line  integration  of  the  transmission 
functions.  Therefore  it  is  necessary  to  make  some  reasona¬ 
ble  approximations  to  compute  transmission  functions. 

The  apparently  random  line  positions  of  the  absorbing 
gases  in  the  absorption  spectra  make  it  possible  to  treat 
absorption  in  a  band  statistical ly.  Goody  (1952,  1964) 
found  that  a  simple  Poisson  distribution  of  line  strengths 
agrees  well  with  experimental  water  vapor  data.  Another 


very  useful  approximation  makes  it  unnecessary  to  integrate 
over  exact  line  shapes,  which  change  with  pressure  and  in 
variations  of  amount  of  absorber.  The  Curtis-Godson  approx¬ 
imation  states  that  transmission  along  a  non -homogeneous 
path  is  very  close  to  that  along  a  path  at  constant  pressure 
with  a  constant  effective  absorber  amount.  As  a  result,  a 
single  line  absorbing  along  a  non-homogeneous  path  can  be 
closely  represented  by  a  line  absorbing  along  a  homogeneous 
path  by  adjusting  the  line  intensity  and  width.  Choice  of 
width  and  intensity  is  made  to  give  correct  absorption  in 
the  weak—  and  strong-line  limits.  Following  the  work  of 
El sasser  (1942)  and  Rodgers  and  Walshaw  (1966)  in  approxi¬ 
mating  the  flux  transmission  function,  a  diffusivity  factor 
of  1.66  has  been  included. 

Mean  transmission  functions  used  for  carbon  dioxide 
were  pre-calculated  line-by-line  by  Dr  Virgil  Kunde  of  the 
NASA  Goddard  Space  Flight  Center  (Krishnamurthy,  1982)  based 
on  an  assumed  standard  temperature  sounding  and  a  uniform 
mixing  ratio  of  CO2.  The  transmission  functions  used  for 
ozone  were  pre-calculated  by  Dr  Noelle  Scott  of  the  Labor a- 
toire  de  Meteorologie,  Paris  (Krishnamurthy,  1982).  The 
multiplicative  property  of  transmission  functions  is  used  in 
case  of  overlap  of  absorption  spectra  of  different  absorbers 
in  a  particular  spectral  interval. 

3.2.3  Clouds 

Clouds  occur  if  the  model  predicts  large-scale  super — 


saturation,  which  consists  of  a  comparison  at  sach  grid 


point  of  ths  saturation  and  actual  mixing  ratios.  When 
clouds  occur  an  iterative  reduction  of  ths  mixing  ratio  (ac¬ 
companied  by  a  rise  in  temperature  due  to  latent  heat  re¬ 
lease)  is  carried  our  until  the  relative  humidity  at  the 
grid  point  equals  100%.  Clouds  are  also  created  when  a  spe¬ 
cific  type  of  comparison  of  pairs  of  layers  shows  the  lower 
layer  to  be  buoyant  with  respect  to  the  higher  layer.  These 


convective  clouds  are  assumed  to  fill  an  entire  horizontal 
grid  box.  However,  since  100%  fractional  cumulus  cloudiness 
is  unrealistic,  radiation  effects  of  convective  clouds  are 
neglected  by  the  model.  This  actually  leads  to  an  improved 
longwave  radiation  distribution  (Moeng  and  Randall,  1982). 
Clouds  are  assumed  to  be  black  bodies,  however,  to  alleviate 
excessive  cooling  due  to  high-level  clouds  in  the  model,  all 
clouds  in  the  top  three  layers  of  the  model  are  considered 
to  be  transparent  with  respect  to  radiation.  Latent  heating 
and  drying  due  to  large-scale  super-saturation  in  these 
three  layers  still  occur.  All  condensate  released  is  assum¬ 
ed  to  precipitate  to  the  next  lower  unsaturated  layer  where 
it  evaporates.  A  cloud  is  created  wherever  condensation 
raises  the  relative  humidity  to  100%.  Each  column  of  layers 
is  tested  from  the  top  downward,  and  condensation  in  layer  9 
precipitates. 


3.2.4  Longwave  Flux 

A  given  volume  of  atmosphere  both  receives  and  emits 


infrared  (longwave)  radiation  and  warms  or  cools  according 
to  the  net  vertical  flux  across  the  volume.  In  computing 
longwave  vertical  fluxes,  the  Wu-Kaplan  radiation  routine 
considers  the  following  absorbers  and  their  associated  ab¬ 
sorption  regions: 


(1)  Water  vapor 


(2)  Water  dimer 


<  0  -  2050  cm-1  ), 

(  660  -  1200  cm-1  ), 


(3)  Carbon  dioxide  (  500  —  800  cm-1  ), 


(4)  Ozone 


<  800  -  1200  cm-1  ). 


Parameters  used  to  compute  longwave  fluxes  at  each  grid 
point  are  pressures,  temperatures,  specific  humidities  and 
cloud  status  at  all  vertical  levels,  as  well  as  latitude, 
number  of  spectral  intervals,  Julian  day,  and  pre-calculated 
transmission  functions.  For  purposes  of  radiation  flux  cal¬ 
culations,  two  stratospheric  layers,  with  upper  boundaries 
at  1  mb  and  5  mb,  are  added  to  the  nine  tropospheric  lay¬ 
ers.  Temperatures  at  the  tops  of  these  two  stratospheric 
levels  come  from  the  climatological  data  of  Newell,  et  al . 
(1974).  In  addition,  each  original  model  layer  is  divided 
into  two  "radiation"  layers  to  enhance  calculation  accura¬ 
cy.  Transmission  functions  and  Planck  functions  are  comput¬ 
ed  at  the  19  radiation  levels.  Use  of  Simpson’s  three-point 
rule,  which  requires  temperature  and  specific  humidity  val¬ 
ues  midway  between  radiation  levels,  requires  that  each  ra¬ 
diation  layer  be  further  divided  into  two  "radiation  sublay- 
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ers“ .  Net  longwave  fluxes  are  computed  at  the  upper  bounda¬ 
ries  of  all  model  layers  and  at  the  lower  boundary  of  the 
bottom  <9th)  layer.  Specific  details  on  numerical  tech¬ 
niques  involved  in  computing  fluxes  are  described  in  Mu 
(1976) . 

3.2.5  Initialization 

The  GLAS  model  was  initialized  with  the  observed  at¬ 
mospheric  conditions  of  15  November,  1978.  Prescribed 
boundary  conditions  included  topography  and  climatological 
values  of  sea  surface  temperatures,  surface  albedoes,  and 
snow/ice  distributions.  Initial  soil  moistures  were  calcu¬ 
lated  from  cl imatological  values  of  surface  relative  humidi¬ 
ties  (Randall,  1982). 


On*  tool  for  objectively  describing  the  GLAS  model's 
skill  in  sieulsting  OLR  is  the  spherical  hareonic  Method  of 
analysis  which  Christ idis  and  Spar  (1981)  call  ua  valuable 
critical  tool  for  coeparing  synoptic  patterns  and  evaluating 
the  output  of  general  circulation  Models." 

4.1  Properties  of  Spherical  Harmonics 

W*  can  express  F,  our  OLR  data  field  (a  function  of 
latitude  and  longitude)  as  a  sum  of  terms  as  follows* 

N 

F(6,4)  *  I  AnPnOl)  4 
n=0 

N— 1  N 

I  I  2  C  An.  mCOS  (m+)  -  Bn.rnSin(m4)  3  Pn.mOD  (1) 
*«1  n*m 

Where  (see  Byerl y ,  1893,  and  Blackmon,  1976)* 

9  *  co-latitude  (90  minus  latitude  in  degrees) 

4  ■  longitude 
M  ■  cosO 

m  -  number  of  waves  around  a  latitude  circle;  order 
of  the  harmonic 
n  ■  degree  of  the  harmonic 
n-m  ■  number  of  zero  points  between  poles 
n  2m  in  all  cases 
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Pn(.M)  *  Ligmdrc  polynomial,  or  zonal  harmonic 
dafinid  as  follows: 


Pn  “  _ I _ S42.<M*-1>"  (Rodrigue's  formula) 

2"n!  dxn 

and, 

Pn.mOi)  m  Associated  Legendre  function  defined  as 
follows. 


Pn.mOl)  *  <1-J4a)  fjm  p„  f  u> 

dwm 


The  expansion  (1)  is  based  on  the  orthogonality 
conditions: 


For  Legendre  polynomials, 

1 

r 

->  PmOOPn,  oodw  *  0  ,  n  A  m 
-1 

and, 

1 

2n+l  r 

->  C  PoOO  3*  dx  -  1 
-1 


Likewise  for  associated  Legend re  functions, 

1 

r 

J  Pn.mOOPK.mOOdX  «  0  ,  n  A  k 
-1 


and, 

1 

H  r 

J  C  Pn.m(H)  3*  dW  -  1 
-1 


with  h  -  I2n+1>  jn-ff? ! 


<n+m) ! 


<A«,  An,*)  and  (Bn,*)  art  the  expansion  coefficients  of  the 


series  where. 


2*  1 
An  -  J  r  r 

J  ->  F(0,*)Pr,  (M)dMd* 
0  -1 


with  J  ■  2n+l  ; 


2X  1 

An, a  »  K  /*  r 

•>  j  F (0,4)Pn. m ()<)cos (m+)d)Jd+ 
O  -1 


2*  1 

Oh  ^  x  ^  f* 

j  j  F<e,+)Pn.m<j»nin<m+)dMd+ 
O  -I 


with  K  -  (2n+l) (n-w) ? 

*  (n+e) ! 


4.2  Application 

In  the  series  (1)  explained  variance  of  the  field  F  in¬ 
creases  as  N  increases.  However,  in  using  a  Modification  of 
coeputer  routines  based  on  code  supplied  by  Dr  Jordan  Alpert 
of  NASA,  (private  comeuni cation)  we  find  that  truncating  at 
N  *  7  (triangular  truncation),  captures  approximately  90X  of 
the  variance  of  the  original  fields.  For  greater  accuracy 
we’ve  chosen  to  truncate  it  N  ■  IS,  enabling  us  to  recon¬ 
struct  OLR  fields  which  explain  at  least  96X  of  the  variance 
of  the  original  fields. 


Sine*  our  two  data  cats  have  different  grid  sizes,  they 


are  not  directly  comparable  for  some  applications.  To  be 
fair  to  the  model  with  its  coarser  grid,  all  fields  of  ob¬ 
served  OLR  have  been  decomposed  into  spherical  harmonics  and 
recomposed  onto  the  same  4  X  S  degree  grid  as  the  model's 
fields.  Specifically,  the  global  average  observed  OLR  field 
for  a  given  month  was  calculated  by  averaging  n  years  of  ob¬ 
served  data  for  that  month,  point  by  point,  to  get  an  avef — 
age  field.  This  average  field  was  then  decomposed  and  re¬ 
constructed  onto  the  coarser  grid.  Calculation  of  the  ob¬ 
served  field's  average  standard  deviation  involved  decompos¬ 
ing  and  reconstructing  each  of  the  observed  months  (9  Janu¬ 
ary's,  for  example)  prior  to  the  appropriate  point  by  point 
summations  over  the  n  years. 

After  representing  both  data  sets  on  the  same  grid 
scale  by  a  sum  of  spherical  harmonic  components,  we  can  look 
at  two  different  aspects  in  our  comparison  of  these  compo¬ 
nents.  The  amplitude  of  a  component  gives  us  a  measure  of 
how  much  of  the  variance  of  the  total  field  is  explained  by 
that  component.  Likewise,  the  agreement  of  the  phases  of 
the  same  components  in  each  data  set  gives  us  correlation 
coefficients  relating  the  two  fields.  The  amplitude  and 
phase  of  the  n,mth  harmonic  are  given  as  follows: 


(AMPn,  m)  *  «  (An,  m)  a  +  (Bn,  m)  3 


< PHASE n, m)  m  tan- 1 (Bn , ./An . m) 
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Calculation  of  correlation  coefficients  and  explained 
variances  by  coeponent  are  described  in  Appendix  I.  Tables 
of  aeplitudes  and  explained  variances  by  wave  number  are 
given  in  Appendices  11  and  III. 


3.  Analysis  and  Results 

In  all  comparisons  presented  hara,  siaulatad  data  rep- 
raaanta  tha  avaraga  by  Month  of  two  annual  cyclas  (i.a,  two 
January’s,  two  February’s,  ate)  producad  by  tha  GLAS  modal. 
Each  month  of  obsarvad  data  shown  is  tha  avaraga  o-f  nina 
yaars  of  that  month,  exempt  -for  February  and  Ha^  which  had 
tan  and  eight  yaars  respectively. 

He  begin  with  subjective  comparisons  o-f  siNwilated  vs 
observed  data  -for  tha  -four  seasons.  Figures  la  and  lb  show 
avaraga  January  OLR,  simulated  and  observed.  Tha  gross  -fea¬ 
tures  of  higher  OLR  values  at  low  latitudes  and  lower  values 
a  high  latitudes  are  reproduced  by  the  model.  The  Miost 
striking  di-f-ference  between  the  two  -fields  is  in  the  magni¬ 
tude  of  OLR.  In  most  places  simulated  OLR  is  considerably 
lower  than  observed.  Also  note  the  stronger  north-south 
gradient  of  simulated  OLR  at  mid-latitudes  of  both  hemi¬ 
spheres.  Finally,  the  model  seems  to  show  less  longitudinal 
structure  in  general  than  observations. 

Figures  2a  and  2b  through  4a  and  4b,  which  represent 
siiNilated  and  observed  fields  for  April,  July,  and  October 
show  similar  disagreements.  The  model  seems  to  underesti¬ 
mate  OLR  the  most  at  high  latitudes  and  to  reproduce  it  most 
closely  in  the  lower  latitudes.  The  only  areas  where  the 
model  overestimates  OLR  are  in  three  areas  of  the  tropiesi 


central  Africa,  Indonesia/Southeast  Asia,  and  along  the 


Observed  January  outgoing  longwave  radiation 


Amazon  River  in  S.  America.  These  areas  are  well  known  re¬ 
gions  o-f  convective  activity  and,  since  the  GLAS  model  does 
not  allow  convective  clouds  to  interact  with  radiation,  it 
Tails  to  appropriately  reduce  the  OLR  here. 

These  qualitative  comparisons  give  us  a  general  idea  o-f 
how  well  the  model  simulates  OLR  but  we  have  another  means 
o-f  comparison  which  is  less  subjective  than  a  simple  visual 
comparison  and  which  gives  us  -further  insight  into  how  and 
why  the  model's  simulation  di-f-fers  -from  observations.  In 
Section  4  we  discussed  how  spherical  harmonics  could  be  used 
to  break  a  -field  down  into  components.  We  now  use  those 
methods  to  compare  some  o-f  the  major  components  o-f  the  two 
data  sets. 

We  know  -from  Figures  la-lb  through  4a-4b  that,  in  gen¬ 
eral,  simulated  OLR  is  too  low.  Figure  S  shows  that  the 
magnitude  o-f  the  difference  in  global  averages  is  about 
20-23  Wm_a  throughout  the  year.  We  also  see  that  the  model 
performs  fairly  well  in  simulating  the  observed  seasonal  cy¬ 
cle  of  average  OLR  with  higher  values  in  the  (northern  hem¬ 
isphere)  summer  and  lower  values  in  winter.  Figure  6  shows 
that,  to  a  lesser  degree,  the  model  also  successfully  simu¬ 
lates  the  seasonal  pattern  of  variation  of  standard  devia¬ 
tion  of  the  geographical  field.  Simulated  standard  devia¬ 
tion  values  range  from  about  8-12  Wm~a  higher  than  observed. 

Appendix  II  shows  the  amplitudes  of  all  spherical  har¬ 
monic  components  (truncated  at  N  -  15)  for  the  four  seasons, 
both  observed  and  simulated.  For  each  month  the  spherical 


Standard  deviation  of.  outgoing  longwave 


v\ 


harmonic  coaponent  corresponding  to  tha  syaaatric  part  of 
wava  n unbar  2  (N  >  2,  H  ■  0)  if  by  -far  tha  doainant  one. 
This  coaponent  is  a  major  contributor  to  tha  spatial  varia¬ 
bility  o-f  tha  OLR  -field  and  is  also  syaaatric  about  the 
equator,  and  provides  a  measure  o-f  the  annual  average.  It 
explains  as  auch  as  2/3  of  the  total  variance  of  the  ob¬ 
served  OLR  field  (see  Appendix  III  for  explained  variances 


by  coaponent).  The  aodel ,  however ,  places  even  more  empha¬ 
sis  on  this  component,  using  it  to  explain  as  much  as  4/5  of 
the  total  variance  of  the  simulated  field.  Figure  7  shows  a 
plot  of  the  simulated  and  observed  N  *  2,  II  ■  0  coaponent 
amplitudes  by  month.  We  see  that  the  amplitude  of  the  simu¬ 
lated  component  is  about  12  Wm_a  too  high  throughout  the 
year.  On  the  other  hand,  the  seasonal  pattern  is  reproduced 
fairly  well. 

The  next  most  doainant  component  corresponds  to  N  *  1, 
M  =  0.  It  is  the  doainant  mode  of  seasonal  variation  being 
antisymmetric  about  the  equator,  and  its  amplitudes  are 
shown  in  Figure  8.  As  in  Figure  7,  the  model  adheres  well 
to  the  observed  seasonal  pattern.  However,  the  magnitudes 
of  simulated  amplitudes  are  much  more  realistic  here,  being 
too  high  by  as  much  as  6  Wm-a  only  in  the  (northern  hemi¬ 
sphere)  summer  and  in  relative  agreement  the  rest  of  the 


year. 


Having  examined  dominant  component  amplitudes,  we  turn 


our  attention  to  the  agreement  between  the  phases  of  the  ob¬ 
served  and  simulated  OLR  components.  Just  as  explained 


variance  Mas  usad  as  a  measure  of  agreement  between  ampli¬ 
tudes,  we  can  use  correlation  coefficients  to  measure  the 
association  between  the  phases  o-f  observed  and  simulated 
-fields  as  a  -function  of  spatial  scales. 

Tables  2-3  show  correlation  coefficients  for  wave  num¬ 
bers  1  through  15  broken  down  into  symmetric  and  asymmetric 
components.  The  symmetric  part,  which  corresponds  to  spher¬ 
ical  harmonic  components  with  N  >  0,  (i.e.  the  zonal ly  aver¬ 
aged  component)  describes  how  well  the  north-south  variabil¬ 
ity  of  the  two  fields  agree  in  phase.  The  asymmetric  parts 
<M  #  O)  measure  the  agreement  of  east-west  variability  (i.e. 
the  departure  from  zonal  average)  between  the  phases  of  the 
observed  and  simulated  fields.  The  tables  are  cumulative  in 
that  the  values  for  N  *  5,  for  example,  represent  a  compari¬ 
son  not  of  wavenumber  5  phases  alone,  but  of  how  well  the 
phases  of  the  two  data  sets  for  waves  1  through  5  agree. 
The  very  high  symmetric  component  correlations  indicate  that 
the  model  does  an  excellent  job  of  simulating  the  phases  of 
the  zonal  OLR  features.  The  phases  of  the  longitudinal  fea¬ 
tures  of  OLR  are  less  accurately  reproduced  by  the  model,  as 
shown  by  the  low  asymmetric  component  correlations.  It  is 
interesting  to  note  that  when  all  15  waves  are  considered, 
the  total  correlation  coefficients  and  the  symmetric  compon¬ 
ent  coefficients  remain  nearly  constant  across  the  seasons, 

while  the  asymmetric  correlation  coefficients  are  considera¬ 
bly  lower  during  April  and  October  than  in  January  and  Ju¬ 
ly.  A  look  at  the  variances  in  Appendix  III  shows  that  a 


Tabla  2.  January  and  April  Correlation  Coafficiants 

JANUARY 


N 

CORRELATION 

SYMMETRIC 

ASYMMETRIC 

COEFFICIENT 

COMPONENT 

COMPONENT 

1 

0.97 

1.00 

0.93 

2 

0.97 

0.99 

0.09 

3 

0.96 

0.99 

0.37 

4 

0.95 

0.98 

0.55 

5 

0.94 

0.97 

0.55 

6 

0.93 

0.97 

0.55 

7 

0.92 

0.97 

0.54 

8 

0.91 

0.96 

0.51 

9 

0.91 

0.96 

0.52 

io 

0.90 

0.96 

0.53 

11 

0.90 

0.96 

0.53 

12 

0.90 

0.96 

0.53 

13 

0.90 

0.96 
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similar  -feature  occurs  in  the  amplitudes.  During  April  and 
October,  both  the  model  and  obser vat ions  show  a  decrease  in 
the  variance  explained  by  the  amplitudes  o-f  asymmetric  com¬ 
ponents. 

Finally,  we  present  the  results  based  on  the  student's 
’ t’  statistic.  Figure  9a  shows  the  difference  between  ob¬ 
served  and  simulated  OLR  for  January.  We  saw  in  Figures  la 
and  lb  that  in  most  areas  simulated  OLR  values  were  lower 
than  observed  values.  Figure  9a  makes  it  much  easier  to  see 
that  the  model  errors  for  January  range  from  60  Wm~a  too  low 
at  some  high  southern  latitudes  and  along  northwestern  N. 
America  to  as  much  as  45  W*-a  too  high  at  the  three  areas  of 
convective  activity  in  the  tropics.  While  these  numbers  are 


useful  in  themselves,  they  become  much  more  so  if 
whether  or  not  they  are  statistically  significant. 


know 


Figure  9b  shows  the  field  of  values  calculated  using  a 
student's  t-statistic  (see  Appendix  I  for  calculations  of 
t -values) .  We  use  these  values  to  test  the  null  hypothesis 
that  the  difference  between  simulated  and  observed  OLR  means 
are  not  significantly  different  from  zero.  The  heavy  solid 
contour,  2.3,  and  the  heavy  dashed  contour,  -2.3,  represent 
the  boundaries  of  areas  of  significance  at  the  .05  level  for 
a  two-tailed  test.  Areas  with  values  between  -*-2.3  and  -2.3 
(i.e.  areas  inside  the  2.3  contour  but  outside  the  -2.3  con¬ 
tour  are  the  only  areas  where  the  means  of  the  two  data  sets 
are  not  significantly  different.  From  Figure  9b,  then,  we 
see  that  nearly  all  values  outside  the  tropics  are 
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statistically  significant  and  values  in  the  three  convective 
centers  within  the  tropics  are  also  statistically  signifi¬ 
cant.  Let  us  now  advance  through  the  seasons,  looking  at 
the  same  types  of  fields. 

Figure  10a  shows  that  the  April  differences  have  moder — 
ated  in  southern  high  latitudes  but  increased  in  northern 
high  latitudes.  Figure  10b  again  shows  values  of  statisti¬ 
cal  significance  nearly  everywhere  outside  the  tropics.  Al¬ 
so,  the  significant  areas  within  the  three  convective  cen¬ 
ters  have  shrunken  slightly.  July’s  difference  field.  Fig¬ 
ure  11a,  shows  a  considerable  decrease  from  April  over 
northern  high  latitudes  in  general  and  over  northwestern  N. 
America  in  particular.  In  addition,  tropical  areas  of  model 
overestimation  have  migrated  to  the  northwest  somewhat  fol¬ 
lowing  the  ITCZ.  The  t-values  for  July  (Figure  lib)  also 
show  a  drift  to  the  northwest  of  statistically  significant 
tropical  convective  areas.  Note  that  the  area  over  central 
S.  America  is  no  longer  significant. 

October’s  Figures  12a  and  12b  show  similar  patterns 
with  a  slight  intensification  of  tropical  areas  of  overesti¬ 
mation  and  once  again  three  tropical  areas  of  statistical 
significance. 
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We  validated  tha  GLAS  6CH’i  OLR  using  observations  with 


high  accuracy  and  global  density.  We  used  conventional  com¬ 
parisons  as  well  as  spherical  harmonic  analysis  and  statis¬ 
tical  procedures  -for  validation.  This  provided  a  measure  o-f 
quantitative,  in  addition  to  merely  qualitative,  assessments 


of  the  model’s  strengths  and  weaknesses  in  simulating  mean 
outgoing  longwave  radiation. 

The  model  produces  low  estimates  o-f  OLR  at  all  lati¬ 
tudes  outside  the  tropics.  This  deficiency  was  previously 
reported  by  Shukla,  et  al .  (1981).  The  basic  problem  seems 
to  lie  in  excessive  model  cloudiness  as  evidenced  by  a  simu¬ 
lated  planetary  albedo  o-f  0.40  (Randall,  1982a)  compared 
with  the  observed  albedo  o-f  about  0.30.  However,  we  can  now 
state  with  more  confidence  the  degree  to  which  the  model  is 
deficient  in  some  areas.  Global  average  OLR  is  too  low  by 
20-23  Wm~*.  This  represents  appr ox i matel y  10%  of  the  ob¬ 
served  global  average  OLR.  Regional  estimates  fall  within  a 
range  from  73  Wm~a  lower  than  observed  to  43  Wm-a  higher 
than  observed.  High-latitude  estimates  in  both  hemispheres 
are  consistently  low  while  OLR  in  the  tropical  convective 
centers  is  routinely  overestimated. 

The  dominant  spherical  harmonic  component,  (2,0),  which 
represents  the  equator-pole  gradient  in  a  global  average 


mom,  shows  us  that  the  OLR  gradient  is  consistently  too 
steep  by  about  12  We-3.  The  component  which  represents  a 
more  seasonal  than  annual  picture,  (1,0),  shows  strong  evi¬ 
dence  that  the. model  successfully  simulates  the  observed 
inter-hemispheric  seasonal  cycle  of  OLR. 

The  model  performs  reasonably  well  in  reproducing  the 
phases  and  amplitudes  of  the  component  waves,  particularly 
the  zonal ly  symmetric  components.  Perhaps  the  lower  values 
of  the  zonal ly  asymmetric  component  correlations  and  simu¬ 
lated  amplitudes  will  provide  some  guidance  as  to  how  best 
to  improve  the  model  dynamics  so  that  east— west  variability 
can  be  better  simulated. 

Statistically,  the  model  fares  rather  poorly  over  most 
of  the  globe  with  the  exception  of  portions  of  the  tropics, 
where  observations  show  variability  of  OLR  to  be  highest. 
The  various  figures  presented  also  show  how  very  poorly  the 
model  performs,  in  a  statistical  sense,  at  high  latitudes. 

The  excessive  tropical  values  of  OLR  generated  by  the 
model  are  surely  due  to  radiative  neglect  of  convective 
clouds.  Perhaps  further  research  along  the  lines  of  this 
study  would  provide  justification  for  re-introducing  code 
into  the  model  to  allow  convective  cloud— radiation  interac¬ 
tion,  at  least  in  portions  of  the  tropics. 

Future  research  could  also  validate  the  year-to— year 

variability  of  the  standard  deviation  of  OLR.  This,  of 

course,  would  require  a  much  larger  sample  of  simulated  data 
than  was  used  in  this  study.  Other  possible  research  areas 


i 


include  validation  of  temperatures,  winds,  stream  functions, 
etc.  We  feel  strongly  that  further  research  efforts  should 
emphasize  objective  and  qualitative  validation  procedures 
similar  to  those  employed  in  this  study. 


I.  Global  Standard  Deviation  (S«> 


N 

S»  «  [  i  (4X/360)  s*(j)  cose,  3*'* 
j  =  l 

Where, 

N  *  number  of  grid  points  between  North  and  South  poles. 

•  j  *  t  90  —  4<j-l)  1  (2X/360)  :  Latitude  converted  -from 

degrees  to  radians. 

M 

sa(j)  -  <1/M>  Z  Z  F<i , j)  -  AVG  3* 
i  =1 

H  *  number  o-f  grid  points  along  a  latitude  circle. 

N 

AVG  *  r  (4X/360)  XBAR<j)cos9j  <  Global  average  OLR. 
j-1 

M 

XBAR(j)  *  (1/M)  z  F(i,j)  t  Average  around  latitude 
circle. 

i=l 

And  F(i,j)  is  our  field  as  a  function  of  i,j  grid  points. 


II.  Monthly  Standard  Deviation  Field  (  S(i,j)  ) 


Y 

S(i,j>  -  C  (1/Y)  Z  <  F(i,j)r,  -  AVGOLR  ( i ,  j )  )*  3*'a 

n«l 

Where, 

Y  *  number  years  of  data 

Y 

■  (1/Y)  Z  F(i,j)r, 
n«l 


AVGOLR (i, j) 
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III.  %  Explained  Variance  <EV) 


EV  *  100<AMPr».*>  */Sa* 


IV.  Correlation  Coe-f -f icients 


A.  Total  Coe-f -f icients  (  R(x,y)  ) 


R(x,y)  *  I  t  Anan  +21  (An  .man,  m+Bn  .Irhn.m)  l/(  s(x)s(y)  ) 
n=0  m*l 


x  and  y  represent  the  observed  and  simulated  -fields, 
respectively. 


For  our  purposes,  N  is  truncated  at  15. 


S(x)  a  C  I  (  Ar,3  +  2  I  <Ar..m*  +  B o.m5*)  >  3  ! 

n=0  m*l 


S(y)  =  C  X  (  an*  +  2  X  (an.m3  +  br..m»)  ) 

n=0  m=l 


B.  Symmetric  Coe-f -f icients  <  R»(x,y>  ) 


Ra(x,y)  *  X  Ana n  /  Cs*(x)  sa(y)l 
n=0 

Where, 


•■(x)  =  C  X  An* 
n«0 


Sa(y)  »  C  X  an* 
n*0 


C.  Asymmetric  Coefficients  (  RA<x,y)  > 


N  n 

RA<X,y>  *  2  X  I  (An, man,  I 
n»l  m= 1 


+  Bn.nbn.m)  /  CSA(x)Sn(y)  ] 


N  n 

C  2  X  X  (An.m2  +  B 
n=l  m=l 


n,*a)  3*'a 


N  n 

iA(y)  =  C  2  X  X  (an, 
n=l  m=*l 


+  bn,ma) 


V.  ’t ’-Values  <t> 


t  =  K  (XBAR  -  YBAR)  /  (n*s».3  +  nvsv3)*'3 


K  =  C  n«nv(n«  +  nv  -  2)  /  (n«  +  nv>  1 


and. 


XBAR  and  YBAR  represent  averages  of  n«  and  nv  years  o-f 
OLR. 

n«  and  nv  are  the  number  o-f  years  o-f  observed  and 
simulated  input. 


*  C  <l/n„)  X  (xi  -  XBAR)3  l1-'3 
i«l 


*  C  < l/nv)  X  <yi  -  YBAR)3  3»'3 
i*l 


Obsarvad  Data,  January 


Observed  Data,  July 


n 

N 

1  13. 

2  20. 

3  7. 

4  7. 

5  2. 

6  2. 

7  6. 

8  4. 

9  3. 

10  O. 

11  1. 

12  O. 

13  1. 

14  O. 
13  O. 


9  10  11  12  13  14 


.22  '■ 

.Ol  0.91 

.40  0. 61  0.43 

.73  0.35  0.58  0.28 

.27  0.19  0.60  0.19  0.33 

.43  0.23  0.27  0.13  0.24  0.49 

.08  0.19  0.34  0.28  0.26  0.34  0.18 

.23  0.17  0.11  0.24  0.26  0.07  0.41  0.32 

.09  0.34  0.35  0.26  0.11  0.17  0.26  0.23  0.06 

.22  0.06  0.19  0.32  0.20  0.22  0.20  0.11  0.29 


Observed  Data,  October 


H 

N 

1  3. 

2  23. 

3  2. 

4  7. 

5  O. 

6  4. 

7  4. 

8  4. 

9  2. 

10  O. 

11  1. 

12  O. 

13  1. 

14  1. 
13  O. 


03  1. 
68  2.' 
48  o.: 
70  i.: 
61  2.: 
Ol  1. 
02  1. 
70  O.! 
20  l.< 
73  i.: 
93  t.< 
37  o.: 
68  l.< 

03  o.: 
17  o.: 


8  9  lO  11  12  13  14 


.39 

.44  2.14 

.01  0.40  0.96 

.63  0.91  0.62  0.61 

.70  1.07  0.29  0.34  0.97 

.63  0.82  0.39  0.58  0.16  0.71 

.73  0.62  0.16  0.49  0.99  0.38  0.44 

.24  0.32  0.26  0.39  0.16  0.23  0.34  0.24 

.44  0.34  0.36  0.40  0.61  0.52  0.21  0.14  0.34 

.19  0.38  0.25  0.33  0.20  0.23  0.30  0.18  0.16  0.06 

.47  0.23  0.46  0.16  0.40  0.43  0.13  0.18  0.22  0.19  0.18 

.61  0.27  0.16  0.15  0.19  0.33  0.26  0.16  0.33  0.02  0.26 


APPEND ix  I»  i  AMPLITUDES 
Simulated  Data,  January 
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0.66 

85.01 

3 
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